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. The asymptotic iteration method is used to find exact and approximate solutions of Schrodinger's 

■ equation for a number of one-dimensional trigonometric potentials (sine- squared, double-cosine, 
I ^ I tangent-squared, and complex cotangent). Analytic and approximate solutions are obtained by first 
^ , using a coordinate transformation to reduce the Schrodinger equation to a second-order differential 

■ equation with an appropriate form. The asymptotic iteration method is also employed indirectly 
to obtain the terms in perturbation expansions, both for the energies and for the corresponding 

a^: eigenfunctions. 
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I. INTRODUCTION 

T— I ■ 
>■ 

J — ■ Analytic solutions of Schrodinger's equation are possible only for few potential models. Although there are a number 

. of approximation methods available in the literature to estimate the solution of Schrodinger's equation, they cannot 

■ usually be used to study both the exact solutions as well as approximate solutions. The asymptotic iteration method 
C^l , was recently introduced and this can indeed find both the exact and numerical solutions of Schrodinger equation for a 

■ very large number of quantum potentials [2] . If the eigenvalue problem is analytically solvable, the aysmptotic iteration 
T-H [ method provides explicit algebraic expressions for the energy eigenvalues and the corresponding wave functions [3|. 

■ Meanwhile, if the eigenvalue problem is not exactly solvable, the asymptotic iteration method provides approximate 
' solutions (for the eigenvalues and eigenfunctions) through its termination condition 

j> ■ 

. , In the present work we use the asymptotic iteration method (AIM) to provide both exact and approximate solutions 

' of Schrodinger's equation for a class of one-dimensional trigonometric potentials. The approach discussed here can 

5—1 ] be easily applied to solve Schrodinger's equation for other potentials, such as the Rosen-Morse potential, the Eckart 

• potential, and the Poschl- Teller potential. 

Another purpose of the present work is to explore perturbation theory based on the asymptotic iteration method to 
obtain a perturbation expansion for the energy eigenvalues P, @] , up to the fourth-order of approximation in some 
cases. The advantage of the perturbation method presented here lies in the direct computation of the coefficients in 
the perturbation series without the necessity of having the eigenfunctions of the unperturbed problem available. 

A variety of applications of the asymptotic iteration method are studied in the present work: first as a method for 
obtaining exact analytic solutions, second as a direct approximation method, thirdly as a method for obtaining the 
terms in perturbation series. The organization of the paper is as follows: in the next section we present a brief 
introduction to the asymptotic iteration method; this is followed by the description of its perturbation version in 
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section III. In section IV, we study the sine-square potential defined by 

We use the standard asymptotic iteration method to obtain accurate numerical approximations to the energy eigen- 
values. We then analyze the eigenvalue problem using our perturbation method, regarding Vi (x) as a perturbation of 
the zero potential. In section V, we study the double-cosine potential defined by 

^ , . ^ _ / ^1 cos(x) -I- V2 cos(2x), < x < 271 , , 

^^^^'~\oo, x<0 and x>2Tr 

by means of the asymptotic iteration method. In section VI, we apply the asymptotic iteration method to obtain 
exact analytical solutions for the tangent potential 

Vsix) ^Vo tan2(-), x € (" ^'TJ ' 
Finally in section VII, we investigate the exact solutions of the complex cotangent potential 

. s f ii^cot(-), < X < air 

V4:ix)=< ° n J (4) 

^ ' I 00, X = (J and X = an ^ ' 

II. A BRIEF INTRODUCTION TO THE ASYMPTOTIC ITERATION METHOD 

A detailed description of the asymptotic iteration method can be found in Refs. 0-0] ■ For completeness, we present 
here a brief outline of the method. Given sufficiently differentiable functions Xo{x) and So{x), the second-order linear 
differential equation 



y"{x) = Xo{x)y'{x) + so{x)y{x) (5) 



has a general solution 



y{x) = exp ^- J adtj ^Ca + J exp (^J {Xq{t) + 2a(T))dr^ d?j 



where Ci and C2 are constants, if for sufficiently large n > 0, 

■^n -^n— 1 



(6) 



(7) 



An A„-l 

where 

An = Ajj_]^ + S„_l + AoA„_i, Sn — s[^_i + SoXn-l, 11=1,2,.... (8) 

Note that the termination condition ([7]) can be written equivalently as 

6„{x) =Xn{x)s„-i{x) - S„{x)Xn-i{x) = 0, (9) 

where n refers to the iteration number. The <5„(x) function has a crucial role for solving the boundary- value problem: 
since the coefficients s„ and A„ depend on the (unknown) eigenvalue E that are usually determined by the roots of 
equation ([9|). If, for Aq and sq and exact E, the termination condition Sn{x) = is satisfied at every x point, then 
the problem is called 'exactly solvable' and in this case the radial quantum number (or the number of nodes in the 
wave function) is equal to the iteration number. If the problem is not analytically solvable, the expression for 6n{x) 
depends on both x and the unknown E where in this case the eigenvalues are obtained through the numerical roots of 
^ at some suitable (chosen) xg-point. This initial value may be determined generally as the maximum value of the 
asymptotic wave function or the minimum value of the potential Q, however, the initial value is easily determined 
if the spacial domain for the given eigenvalue problem is bounded. Indeed, for this class of problem where x G [a, b], 
the average point xq = {a + b)/2 will be sufficient to initiate the iteration process. 
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III. A PERTURBATION THEORY VIA ASYMPTOTIC ITERATION METHOD 

In Ref. the present authors gave a new perturbation theory based on the asymptotic iteration method @. In 
order to apply the method for a class of trigonometric potentials regarded as perturbations of the zero potential, we 
outline our perturbation approach. Assume the potential is given as follows 

V{x) ^Vo{x)+ fiViix) (10) 

where Vo{x) is a solvable potential, Vi{x) is the perturbation potential, and is a perturbation expansion parameter. 
The one-dimensional Schrodinger equation reads 

-^ + {Voix)+iiVi{x))^ix)^E-i'ix). (11) 

Substitute 'i!{x) = y{x)f{x) where y{x) is an asymptotic solution satisfying the boundary conditions for the eigenvalue 
problem ([TT]). Substitution yields the following differential equation for f{x) 

f"{x) = Ao(a;, fi)f'{x) + so{x, fi)f{x), (12) 

which is in a form suitable for the asymptotic iteration method. Using (jl2p . we can now compute the sequences 
Xn{x,n) and s„(a;,/i). The termination condition Sn{x,^), as given by (fT2|) . can be written as 

5n{x;H,E) = \n{x] ^1, E)Sn~l{x\ IJL, E) - Sn{x\ H, E)\n-l{x\ E) = 0. (13) 

Expanding Sn{x; fi, E) about /i = 0, we obtain the following infinite series 

CO 

E)^J2^^^^kix■,0,E)■^i'' (14) 

fc=0 

that resembles Taylor's expansion with 

1 d^5n{x;iJL,E) _ 1 ^^[\n{x]^J.,E)Sn-l{x■,|J^,E) - Sn{x; IJL,E)Xn-l{x\ jJL.E)] 



5n,k{^,Q,E) 



, fc = 0,l,2,.... 

According to the standard theory of the asymptotic iteration method, the eigenvalues are obtained by the solution of 
the equation (5„(x; fi, E) — for arbitrary fi and thus the perturbation coefficients satisfy 



Sn,kix;0,E)=0. (16) 

It is clear from the standard theory of the asymptotic iteration method that the root of the equation Sn_o{x; 0, E) — 
gives the eigenvalue of the non-perturbative Hamiltonian, say, En_o that may depend on n (the number of iterations 
used by the asymptotic iteration method). Meanwhhile Sn,i{x;0, E) = gives the first-order correction to the 
eigenvalue En^, say En^i, 6n,2{x] ij.q,E) — gives the second-order correction to the eigenvalue i?„,o, say i?ra,2, and 
so on. This procedure allows us to find the coefficients of the eigenvalue expansion 

oo 

^ = (17) 

fc=0 

where (again) En^k are computed by solving Sn^k{x,0, E) — 0. It is an attractive feature of the asymptotic iteration 
method that we can also find the eigenfunctions by means of equation ^ where it is sufficient to compute 



/(x)=exp|^-y a{t,fi)dtj, a{x, fi) = Sn{x, fj.)/ X„{x, fi) (18) 
Again if we expand a{x,ii) about /i = 0, we obtain the following series 

oo 

a{x, A) = ^ ak{x, 0)^*", (19) 



fc=0 
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where ak{x, 0), fc 0, 1, 2, . . . is given by 



akix,0) 



1 d''a{x,X) 



The correction of the wave function is then given by 



(20) 



-E^y afc(T,0)dTj = =expf-/y «fe(T,0)drj, fc = 0,l,2,.... (21) 

Thus we have a perturbation method (based on the asymptotic iteration method) that could in principle provide a 
complete solution to the given eigenvalue problem. 

IV. SINE-SQUARED POTENTIAL 

In this section we present the solution of the one-dimensional Schrodinger equation for the trigonometric sine-square 
potential H]) 



^{x)=E^{x\ vl;(!^) =v|/(_!^) =0, a^O. (22) 



The behavior of the wave function at the boundaries suggests a solution of the form 

^{x) - cos (^) !{x). (23) 
Substituting for "^{x) in equation ^F^ . we obtain 

rW^^.an(;)fW+(Jj-£ + K.si„»(i))/W 

A further change of the independent variable x — a arcsin(?/) yields 



(24) 



where u = 1 + /i — a^E and = Voa^. We now find solutions of this differential equation by using two different 
approaches: (1) the standard asymptotic iteration method discussed in section II; (2) the perturbation method 
developed in section III. 

A. The standard asymptotic iteration method 

We may initiate the iteration method with 

The sequences s„ and A„, n = 1,2, .. . can then be evaluated using These quantities, for fixed a and Vb, are 
functions of the parameter E and the variable y. For a suitable initial value oi y — < 1, the eigenvalues E are 
computed using the roots of the termination condition We select the initial value j/o = that locate potential 

minimum as given by ([T|) . In Table 1 , we report our numerical calculation of En for different values of the parameter 
where n represent the node of the wave function solution, as given using the asymptotic iteration method. 
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TABLE I: The computation of a^E where /j, — Voo^ using the standard asymptotic iteration method. The subscript refer to 
the number of iteration. 



n 


a^Eifi = 0) 


a'^Eifi = 0.1) 


a^E{fi = 0.5) 


a^Eifi = 1) 


a^E{n = 5) 


a'^Eifi = 10) 





1 


1.024 922 118 9io 


1.123 077 224 814 


1.242 428 826 O14 


2.082 985 293 220 


2.923 668 494 222 


1 


4 


4.049 947 916 8ii 


4.248 698 005 O14 


4.494 793 078 615 


6.370 661 125 O21 


8.492 474 366 723 


2 


9 


9.050 038 818 612 


9.250 946 208 9i6 


9.503 664 867 Oie 


11.569 339 156 922 


14.185 709 970 I24 


3 


16 


16.050 020 833 2i3 


16.250 520 743 817 


16.502 081 901 Oi7 


18.551 201 398 423 


21.194 837 346 927 


4 


25 


25.050 013 020 814 


25.250 325 524 he 


25.501 302 132 2i8 


27.532 566 336 24 1 


30.130 124 200 228 


5 


36 


36.050 008 928 615 


36.250 223 215 3i7 


36.500 892 873 821 


38.522 331 587 4 


41.089 436 879 727 



B. Perturbation theory via the asymptotic iteration method 

The exact values of a^E in the case = as shown in Table 1 do not come as a surprise for we have noted the 
problem is equivalent to solving Schrodinger equation for the particle in rigid box Vq — 0; namely for the potential 



f X e (-— —) 

V(x] = < ■^'^V 2'2/ 

^ ' [ C50 l^l > ^. 

The differential equation reads for /i = 

(l-y')/"(2/)-3y/'(y)-c./(y) = 
can immediately solve with polynomial solutions are given for 

to = -n{n + 2)^ a^E ={n + if 

by 

fniv) - 2Fi{-n, « + 2; I ^ /„(x) = 2Fi[-n, n + 2; | i - ^sin(^)). 



(26) 

(27) 
(28) 

(29) 



1, Energy Expansion 



Suppose that the a^E can be written as follows 



k=0 



(30) 



In order to get a zeroth-order approximation to a^E using the iteration method discussed in section 3, we have for 

<5o(x,0) = (31) 

that exact values 

i^n,o = (f^+1)', n = 0,l,2,... (32) 
as expected using the exact solutions mentioned earlier. Thus for the ground state n = 0, we have 

o^Eq = 1 + i/Qj • + t^o,2 • /i^ + . . . (33) 
In order to find the first order approximation I'o,!, we take o?Eq = 1 + 1^0,1/^ ^^nd solve the termination equation 



^i(y,0) 



9/i 







(34) 



/j=0 



where S{y^ fi) is now computed using 



Ao 



3y 



„ _ -i^Q,iA'+My 
So - fZip — 



Straightforward computation using the sequences A„ and Sn, we obtain for the second-iteration the following expression 
for S{y,fi) 



Thus 



Siiy,0) 



(y2 _ 1)3 



(35) 



as the first correction for the eigenvalue. A similar computation using o^Eq — l + ^/i+!^o,2M^ shows that z^o,2 = ^1/128, 
and, using a'^Ea = 1 + (1/4)^ - (l/128)/i2 + l'o,3A^^ that z^o.s = 1/4096, and for q'^Eq = 1 + (1/4)^ - (l/128)/i2 + 
(l/4096)/i'^ + vo^AiJL^ that j^o,4 — —1/393216, where in each case we update the expressions for Ao and sq to initiate 
the computation of (5(y,/z). Thus the ground-state energy eigenvalue expansion takes the form 



1 



1 



a?' En = 1 H — II t 

" 4' 128^^ 4096 



1 



393216 



(36) 



for sufficiently small as usual for standard perturbation theory. For the first-excited state, we obtain the eigenvalue 
expansion 



a^Ei = 4 + i/i ^—i? 

2' 192' 



(37) 



and for arbitrary n > 2, 



a'E^ = (n + 1)2 + ^ + 



2 32n (n + 2) 



(38) 



TABLE II: The computation of a^E for given /i = Vofl^ by perturbation theory via the asymptotic iteration method. 



n 


a^Eifj, = 0.1) 


a'^E{fi = 0.5) 


a'^Eifi = 1) 





1.024 922 118 9 


1.123 077 233 6 


1.242 429 097 5 


1 


4.049 947 916 7 


4.248 697 916 7 


4.494 791 666 7 


2 


9.050 039 062 5 


9.250 976 562 5 


9.503 906 250 



2. Wave function expansion 



In order to calculate the perturbation expansion for the wave function, we note first that the zeroth-order solution of 
the Eq.([22|). using equation ([29]), is 



fn{y) ^ 2Fi{-n,n + 2;^;^-^), n = 0,l,2,. 



(39) 



Unlike the standard perturbation theory, the first-order corrected wave function using our approach can be computed 
without the use of the zeroth-order solution. Indeed, from (|20l) we have for the first correction, using a^E = 1 + t'ci/x. 



a{y,n) 



fiy{-2 + biyoA - Sy^) 



3 -I- 12y2 + - jy^^^t + fii^o jy'^ - 



(40) 
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Thus 



M=0 



and 



^ 9a(r, ^) 



9// 



y(-2 + 5t/o,i - 3y^) 
3 + 12y2 



rf-=^ + ^ln(l+V)~^ln(l + V). 



However, since z/o,i = \ we have 



Finally, the first-order approximation for the wave function is 

h{y) = exp(-^y2) ^ ^ e-t-"'(^). 



(41) 



(42) 



(43) 



(44) 



where is the first-order approximation of the ground-state wave function solution of (j23p . Similarly, we can 

compute the first-order approximation for n — \ where we have 



A(y) = exp(-^y2). 



the first-order approximation for n = 2, we have 



/i(y)=exp 



y 



16 ' 64(4y2 - I] 

Finally, for the first-order approximation in the case of n = 3, we have 



fiiv) = cxp - 



y 



1 



20 60(22/2 _ 1) 



(45) 



(46) 



(47) 



C. Equivalent potentials 

We may note using the trigonometric identities 

sin^ (61) = 1 -cos^ (61) and sin^(6') = ^ (1 - cos(26')) , sin2(0) = cos2(e') - cos(26l) (48) 

that the eigenvalue solutions of the following schrodinger equations for trigonometric potentials confined between two 
rigid walls of infinite potential energy 

*.(^^Ty)=0' (49) 
*«(^^T^)-0, (50) 

*c(-^Tf)=0, (51) 

are related to the eigenvalues of the sine-square potential as given by equation l|22p . These relations are easily derived 
and are given explicitly by 

E^ = E-Vo, E^^E-^Vo, Eq = E. (52) 



V.cos>(i) 



d' 1 (2x 



Vb cos^ f — ^ — Vq cos 



*ax) 

vfclx) 



£;^*^(x), 

E^^^{x), 



8 



V. DOUBLE-COSINE POTENTIAL 

In this section, we consider a single particle in the double cosine potential [loj . 



Vi cos(a;) + V2 cos(2x), < a; < 27r 
00, X < and x > 2Tr 



and obeys the Schrodinger equation 



+ Vi cos{x) + Vi cos(2x) 



'<\)(x)^E^(x), (0) = (27r) = 0. 



(53) 



(54) 



The maximum of the potential (|53p occurs within the given interval between x — Q and x = 21: and has the value of 
Kriaa; = ^1 + V2, whilc the minimum occurs a.t x = tt with the value of Vmin = V2 — Vi. Clearly, if Vi = and V2 > 0, 
the problem corresponds to Mathieu's equation jl^l- Assume the general solution of the differential equation ([5^ 
that satisfies the boundary conditions takes the form = sin(a;/2) f{x). Then direct substitution, using equation 
I, yields the second order differential equation for f(x) as 



f"{x) = - cot (I) fix) + {\ + V, cos{x) + 2V2 cosHx) -¥2- E)f(x) 



Assuming, further, that 



f{x) = exp 



cos{x) g{x), 



we obtain 



g"{x) = (cot (I) + V, sin (x)) g'ix) + (^i - }^ - Yl + (^Yl + 2V2) cos^lx) - V2 - g{x) 



(55) 
(56) 
(57) 



where we have used the trigonometric identity cos(a;) = — 1 + 2cos^(a;/2). Finally with the substitution y = cos(a:/2) 
we obtain 



'(y) = - 



V 



Clearly we have the exact solution 



- ^{yt + W2)y' g{y) 



9oix) = 1, 



subject to V2 = -^Vx 
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that is 



1 Vi 



4 2 

For excited states we have applied the asymptotic iteration method that we used with 



4^12/ 



(58) 

(59) 
(60) 

(61) 



So 



l-2\/i+4\/2-4g 



and a starting value of = corresponding to the minimum of the double-cosine potential at Xq = tt. In Table 
ml we report the first few excited states along with the iteration number used by the asymptotic iteration method. 
All the eigenvalues are accurate for the given number of decimals shown and can easily be obtained with any desired 
precision. 

This work can be easily extended to study the Schrodinger equation with the periodic potential 



dx'^ 



+ vi cos 



V2 COS 



^x) = £:f„^^„^jV(a^), 0<x< 2fc7r, fc = 1, 2, . . . , ^ (0) = V (2M = 0, (62) 



with rigid walls at the boundaries. A simple scaling argument x ^ ax with a = k shows that 



dx'' 



+ k^vi cos(a;) + k'^V2 cos(2x) 



ip^k^Eil), < a; < 27r, fc = l,2,..., * (0) = * (27r) = 0, 



so that we can easily compare with our earlier work concerning equation (|54p. where in this case 

E{k\x,k\2)=k^£l,,,,y 



(63) 



(64) 
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TABLE III: The eigenenergies of Schrodinger's equation with the double cosine potential V{x) = cos(a;) — |cos(2a;), where 
yo = 0. 



n 


En 


n 


En 





- 0.375 000 000 OQQ3 


5 


6.270 895 381 O5I52 


1 


0.625 000 000 OOO50 


6 


9.014 386 533 89O55 


2 


1.258 305 195 0635i 


7 


12.260 510 491 I8856 


3 


2.311 907 935 56452 


8 


16.008 020 136 79557 


4 


4.031 642 235 I9253 


9 


20.256 322 480 09258 



VI. TANGENT-SQUARED POTENTIAL 

Recently, Ta§eli [11] studied the squared tangent potential V{x) = v{v — l)tan^(a;), x £ (— 7r/2,7r/2) also known as 
the symmetric Poschl- Teller potential It was noted that the energy eigenvalues, but not the eigenfunctions, 

of the Hamiltonian with the squared tangent potential on the symmetric interval (— 7r/2,7r/2) are precisely the same 
as those of the Hamiltonian with the squared cotangent potential V{x) — v{v — 1) cot^ x on the asymmetric interval 
(0, tt) considered by Marmorino in fl^ (see also [3). This observation can be easily verified using the identity 

tan(a; + ^) = ^ cot(x). (65) 

In this section, we will used the asymptotic iteration method to study the Schrodinger equation with the confining 
potential 



V{x)=Votan^[-j, a; e — , — j , Vo > 



contained within an infinite square well with sides at x — ± TTa/2 and Vq gives an indication of how rapidly the 
potential increases within the well. In this case, the time-independent Schrodinger equation reads 



(f ■) /x\ 



^ix)^E^{x), VI, _ = 0. (66) 



dx^ 

We, however, take a slightly different approach to solve (|66l) by letting the wave function assume the form 

*(x) =cos(^) ./(x), (67) 



that satisfies the boundary conditions of (|66l) . Again it is straightforward to verify that the substitution of (|67l) into 
(pS)) . yields a differential equation for f{x) as 



fix) = - tan(-)/'(a;) + { - - E + Vo tan2(-) f{x). (68) 
a a " 

A further change of variable x — a arcsin(j/) implies 



f{y), (69) 



and, with the substitution 

/(2/) = (l-2/'r5(y), (70) 

we finally obtain for g{y) 

(4a + 3) (^ + 2a) 
Ti 2r 9 (y + T^ 2V 

(1-y) (1-2/) 



ny)-'^^^^9y)+'^^9iy), y€(-i,i). (71) 



where we set lu = \ — a?E, ji ~ o^Vq and a = (\/ /i + 1 — l) /2. This differential equation is an example of an exactly 
solvable problem solved by means of the asymptotic iteration method. Indeed, with 

(4a + 3)y uj + 2a 
Ao(y) = — ; — and SQ{y) = 



the termination condition (|9l), 5n = yields 

(5i = =^ (w + 2a){w + 6a + 3)=0=^w = 
52=0=^ {w + 2a){w + 6a + 3)(w + 10a + 8] 



—2a, w = —6a — 3 
= w — —2a, w 



-6a — 3,w = —10a 



and in general 

n 

(5„ = =^ (w + + 4a + 2) + 2a) = =4> «;„ = -n{n + 4a + 2) - 2a, 

Thus the energy eigenvalues are given by 

a^En =n^ + {2n + l)(2a + 1), n = 0, 1, 2, . . . 

For the analytic solution of gn{y), we note by equation (6) that 

n = 0=^ .9o(y) = 1 

n = l=^ gi{y) = y 

n = 2=».g2(y)=4(a + l)y2-l 

" = 3 =^ 53^) - y{2{2a + 3)2/3 _ 3) 

n = 4 =^ gi{y) = 8(a + 2)(2a + 3)/ - 6(6 + 4a)y2 + 3 



n = 0,1,2,. 



These polynomial solutions are easily generated by 

9n{y) = 2Fi{~n, n + Aa + 2; 2a + | ^-^), n = 0, 1, 2, . . . 
up to a multiplicative constant. Here, the Gauss hypergeometric function 2^'i(a, /3; 7; x) is defined by [l] 

a P 
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_ {a)kil3)kx'' _ aP X a(a + l)/3(;3 + 1) a;^ 
"^'"^ h)k "fcl^^^Tl!^ 7(7+1) 2! 



fc=0 



and (A)fe denotes the Pochhammer symbol defined, in terms of Gamma functions, by 

r(A + fc) _ J 1 if (fc 0; A e C\{0}) 



(A) 



k • = 



r(A) [ A(A + l)(A + 2)...(A + fc-l) if (fceN;AeC). 
Thus the exact analytic solutions of Schrodinger's equation (|66p is 

*„(a;) = An cos2"+i (f ) 2i^i(-n, n + 4a + 2; 2a + |, i - i sin (f)), n = 0, 1, 2, . . . 

^ a2£:„ = n2 + (2n+l)(2a + l), a = i {y/a^Vo + 1 - l) . 

where An is a normalization factor. Clearly, for Vq = we obtain the exact solutions for the rigid box ([29 
the trigonometric identity 

tan^(x) = scc^(a::) - 1 
we can easily obtain the energy eigenvalues for the Schrodinger equation 

d2 



+ Vo sec" (^) ^-(2;) = £,ec^(,/,)*(x), * (-^) = * (^) - 0. 



where 



^SGC^ (a; / a) -^tan' 



^(.) + K) = ^ (n' + (2n + l)Va'Fo + l + aVo) 
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Using a simple scaling argument x — ^ ax, we can write (|77p as 



+ Uq sec (x) 



where Uq = o^Vq and 2; = a'^£scc^(x/a)- Using the identity 



sec(x + 2'-' ~ ~ csc(a;) 



we can write (l79l) as 



ip^€ip, V(0) = V'W = 0, e = ?i2 + (2n + l)VC^+T+C/o. 



(79) 



(80) 



(81) 



VII. COMPLEX COTANGENT POTENTIAL WELL 

In this section we discuss a quasi-exact trigonometric potential model with complex coefficients 0], 

N I iwcot(-), 0<x<aTT, V > 0,i = 
00, X — and tto 



(82) 



We show that the asymptotic iteration method can provides the necessary conditions under which the one dimensional 
Schrodinger equation 



d'^ / x\ 

-—5- + zi; cot ( - ipix) — Eip(x), ip (x ^ 0) = ^Jj (x ^ an) = 
dx'' V a / J 

has exact (analytic) solutions. To this end we assume that the wave function 'il'{x) takes the form 

V'(x) = sin (^^^ f{x), 

which yields, up on substituting into equation (j83p . the following differential equation for f{x) 

fix) = cot (-) fix) +(^-E + Vo cot (-)) fix) 

wherein Vq = the change of variable y — cot(2:/a) allows us to write (j85p as 

'1 - a^E + a^Voy" 



f"{y) 



f{y) 



(l + y2)2 
and finally, using the substitution 

/(y) = e-^ arccot(,)(^ ^ 2/')"/'5(y), -c^ < 2/ < 00 
where a and /3 are constants to be determine shortly, we obtain 

2iay + (3) a(Q - l)y^ + (2/3( a - 1) - Uoa2)j/ + a + £;a2 + /J^ - 1 

9 (y) = — 1^^2 g (y) 

To solving this equation we note first that 



(1+2/2)2 



iv) 



(83) 



(84) 



(85) 



(86) 



(87) 



2(a-l)' 



a ^ 1 



which yields 



9"{y) 



2{ay + P) , / "'^"(a-i)'"'^ , 

—5 (y) - a(a - 1) I ), , I .9(y) 



1 + r 



(l + y2) 



(89) 



(90) 
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Thus with 

^2 771 /„. -1 \2 q2 



a'E = {a-iy - (91) 
we obtain a more compact form of the second order Unear differential equation for g(x) as 



„ 2{ay + l3) , a{a - 1) 

1 + y 1 + y 



The asymptotic iteration method can then be apphed with 



Ao 



Q(a — 1) 



(93) 



to yield, by means of the termination condition ^ , the following condition for an exact solution 

(5„ = =^ ^— — — ^ TT(a + fc)2 = ^ a = -n, n=l,2,... for a ^ 1. 

a + n 

k=l 

The (analytic) polynomial solutions are given, using ([6]), explicitly as 

71 = =^ go{y) = 1 

71 = 1 =^ .91 (y) = -2?/ + 2/3 

11 = 2 =^ .g2(y) = - 6/3y + 2/?^ _ ^ 

71 = 3 =^ 33 (y) = -3i/ + 9/?y2 - (6/?^ - 3)y + - 2/3 

n = A=^gi{y)^ 15y^ - 60f3y^ + 30(2/3^ - l)y^ - 20P{P^ - 2)x + 2^^ - 10/3^ + 3 

The recurrence relation to generate these polynomial solution is given by 

_ / 2(271+1 + 2a)(7i^ + 1 + a)y 2/3(271 + 1 + 2a){a - 1)\ 4(7i + l)(7i^ + 1 + a)(7i2 + 2an + a'^ + /J^) 

^"+' " ^+2^^ ^ (n + 2a)(7i + a) j (71 + 2a)(7i + a) ^" 

(94) 

initiated with 

go = l, gi=2ay + 2/3 (95) 
In terms of the hypergeometric functions, the polynomial solutions can be written, are 

gn{y) = C„(-2i)" {a + i/3)„ 2F1 ( -n, 71 - 1 + 2a; a + 7/3; i - ^ j , 71 = 0, 1, 2, . . . (96) 



The normalization condition C„ can be computed using 

, 1 (—a + /J7 + 1)1 (—a — pi + 1) [2a)n (za + zri— 1) 

which is valid for a < i — 71, ri = 0, 1, 2, . . . , where Sum — for n ^ m and Snm = 1 for n ~ m. 



(97) 



VIII. CONCLUSION 



We have applied AIM to solve Schrodinger's equation with a number of confining trigonometric potentials. When 
a termination condition is satisfied, the method yields exact analytical solutions. In other cases approximations are 
found either by a direct application of AIM or by using the method to find exact expressions for the terms in a 
perturbation expansion. In particular, our results lead to a new exactly solvable problem with a complex cotangent 
potential which is infinite outside the interval (0,7ra). 

The asymptotic iteration method can be applied to a wide variety of one-dimensional or spherically symmetric quantum 
potentials defined on infinite or semi-infinite domains. In the present work, the method is applied to a number of 
quantum (trigonometric) potentials defined on finite domains [a, b], for which the method proves to be extremely 
effective for finding both analytical and approximate solutions. It is also clear from the present work that the method 
is applicable to periodic potentials and it is our intention to treat such problems in future work. 
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